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Abstract 

Since the energy momentum tensor of a magnetic field always contains a spin-2 component 
in its anisotropic stress, stochastic primordial magnetic field (PMF) in the early universe must 
generate stochastic gravitational wave (GW) background. This process will greatly affect the 
relic gravitational wave (RGW), which is one of major scientific goals of the laser interferometer 
GW detections. Recently, the fifth science (S5) run of laser interferometer gravitational-wave 
observatory (LIGO) gave a latest upper limit 0,gw < 6.9 x 10^^ on the RGW background. 
Utilizing this upper limit, we derive new PMF Limits: for a scale of galactic cluster A = 1 
Mpc, the amplitude of PMF, that produced by the electroweak phase transition (EPT), has to 
be weaker than i?A < 4 x 10^^ Gauss; for a scale of supercluster A = 100 Mpc, the amplitude 
of PMF has to be weaker than Bx < 9 x 10~^^ Gauss. In this manner, GW observation has 
potential to make interesting contributions to the study of primordial magnetic field. 

PACS numbers: 95.85.Sz, 04.30.Tv, 04.80. Nn 
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1. Introduction 



Although the magnetic fields in galaxies and galactic clusters have been observed for many years [1], 
the origin of these magnetic fields still remains a mystery. One of the most promising candidate origins is 
the primordial magnetic field (PMF) produced in the early universe [2]. Such a seed field can be produced 
by the inflation [3, 4, 5], or by the electroweak phase transition (EPT) [6, 7, 8]. In the past, PMF has been 
constrained mainly by using their various effects on cosmic microwave background (CMB) anisotropies 
and polarization [9, 10, 11, 12, 13, 14]. Based on the Faraday rotation of CMB linear polarization caused 
by PMF, Kahniashvili et al [15] utilized the 5-year WMAP (WMAP5) S-mode polarization limit to give an 
upper limit on the amplitude of PMF: Bx < Q x 10^^ Gauss on a scale A = 1 Mpc. Besides, Kristiansen 
and Ferreira [16] forecasted the constraints on PMF from future CMB polarization experiments. On the 
other hand, a stochastic PMF itself can generate the relic gravitational wave (RGW) background [17, 18). 
This is because the energy momentum tensor of a magnetic field always contains a spin-2 component in its 
anisotropic stress. Utilizing the nucleosynthesis limit on gravitational wave (GW), Caprini and Durrer [19] 
gave extraordinarily strong constraints on the amplitude of PMF: for PMF produced by EPT, Bx < 10~^^ 
Gauss; and for PMF produced by inflation, Bx < 10~^^ Gauss. In addition, Giovannini and Kunze [20] gave 
a comprehensive treatment of the PMF's effects on CMB maps, and obtained a constraint Bx<5x 
Gauss by using the WMAP5 B-mode polarization limit. 

Laser interferometer GW detection is one of the most important approaches to detect the RGW. Using 
LIGO S4 data and cross-correlation technique [21], LIGO scientific collaboration [22] gave an upper limit 
^GW < 6.5 X 10~^ on the total energy density of RGW. In addition, LIGO scientific collaboration and 
ALLEGRO collaboration [23] also gave an upper limit Clgyj{u) < 1.02 on the stochastic GW background. 
Especially, in a latest work, using LIGO S5 data, LIGO scientific collaboration and Virgo collaboration [24] 
gave a latest upper limit ^gw < 6.9 x 10^^, which improves LIGO S4 result for one order of magnitude. 
So in this work, we shall give the constraints on PMF by using the latest LIGO S5 result. For the PMF, 
we shall focus on that produced by EPT. As seen below, it yields a level spectrum of RGW in the range 
(10^^ — 10^) Hz, which covers the bands of operation of LIGO. PMF can also be produced by inflation. 
However, since the PMF produced by inflation affects the RGW only in very high frequency range v > 10^ 
Hz that is out of the bands of operation of LIGO, we will not discuss this class of PMF. Our work differs 
from previous papers in the following two aspects: First, in previous papers [19, 20], PMF is directly 
treated as the source of RGW. However, RGW has already existed before the EPT. That is to say, if one 
treated the PMF produced by EPT as the source of RGW, a lot of important information about ultra-early 
universe, such as inflation and reheating, will be lost. Therefore, for a complete treatment, the PMF 
produced by EPT should only be viewed as a subsequent factor that affect the RGW, which is adopted 
in our work. Second, in previous papers, only the nucleosynthesis limit [19] and the WMAP5 B-mode 
polarization limit [20] have been used to constrain the RGW and to give the constraint on PMF. Since the 
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laser interferometer GW detection is also a crucial approach to detect RGW, it would be very interesting 
to give a constraint on PMF by using the data of current GW detection. To our present knowledge, this 
issue has not been investigated before. So in this work, we shall give new limits on the amplitude of PMF 
by using the latest LIGO S5 data. 

The organization of this paper is as follows. In section 2, from the inflationary stage up to the 
accelerating stage, the scale factor a(r) is specified by the continuity conditions for the subsequential 
stages of expansion. In section 3, the analytical solution of RGW is determined with the coefficients being 
fixed by continuity condition joining two consecutive expansion stages. In particular, the important impact 
of PMF produced by EPT is included into our calculation. In section 4, we calculate the spectrum h{u,To) 
and the one-sided spectral density S'/i(z^) of RGW, and then give the theoretical predictions of RGW, which 
can be directly compared with the latest LIGO data. In section 5, by comparing the theoretical predictions 
of RGW with LIGO S5 data, we give the constraints on the amplitude of PMF. Several technical derivations 
are deferred to two appendices. It should be noted that, our conventions are different from that of Ref. 
[19]: the scale factor a(r) has the unit of cm, while the conformal time r and the co-moving distance x 
do not have units. Besides, Greek indices run from to 3, Latin indices run from 1 to 3, the subscript 
"0" always indicates the present value of the corresponding quantity, and the unit with c = h = ks = ^ 
is used. 



2. Expansion History of The Universe 

From the inflationary stage up to the currently accelerating stage, the expansion of the Universe can 
be described by the spatially flat Robertson- Walker spacetime with a metric 

ds^ = a2(T)[-dr2 + Sijdx'dx^]. (1) 

The scale factor a(r) for the successive stages can be approximately described by the following forms [25]: 
The inflationary stage: 

a(r) = Zi|r|^+^, -oo < r < n, (2) 

where 1 + (3 < 0, and ri < 0. This generic form of scale factor is a simple modelling of inflationary 
expansion, and the index /3 is a parameter. If the inflationary expansion is driven by a scalar field, then 
the index (3 is related to the so-called slow-roll parameters, 77 and e [26], as (3 = —2 + {rj — 3e) . in this 
case one usually has (3 < —2. In addition, the WMAP5 data combined with Baryon Acoustic Oscillations 
and Type la supernovae give = 0.960lo:oi3 (95% CL) [27]. Base on the relation = 2/3-1-5 [28, 29], 
the inflation index P = -2.02 for = 0.960. 
The reheating stage: 

a(r) = a^|r-rp|^+^% n < r < r^, (3) 
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where Tg is the beginning of radiation era and Tp < ti. As usual, the model parameter will be taken as 
Ps = -0.3 [30, 31]. 

The radiation-dominant stage: 

a{T) = ae{T - Te), Tg < T < T2. (4) 

This is the stage during which the EPT produces PMF. We use Tew to denote the starting time of EPT: 
Tg < Tew < T2- The Corresponding energy scale is T ~ 100 GeV. As seen below, the wave equation of 
RGW is still homogeneous for r < Tew, but becomes inhomogeneous for Tew <t <T2- 
The matter-dominant stage: 

a{r) = am{T - Tmf , T2<T < Te, (5) 

where < T2 and te is the beginning time of the acceleration era. 
The accelerating stage: 

a{T) = Ih\t - Ta\~'^ , TE<T<To, (6) 

where tq is the present time and tq < Ta- The index 7 in Eq.(6) depends on the present fractional dark 
energy By fitting with the numerical solution of the Friedmann equation [29, 30, 31], 

/a'\2 SttG, , 

where a' = da/dT, one can take 7 ~ 1.044 for $7a = 0.75. The redshift of the start of this stage depends 
on the specific models of the dark energy. For instance, in the KCDM model with JIa = 0.72 and 
rim = 0.28, it starts at 2; ~ 0.37; and in the 3-loop Yang-Mills condensate dark energy model [32], it 
starts at z ~ 0.5. 

In the above specifications of a(r), there are five instances of time, ri, Tg, T2, te, and tq, which separate 
the different stages. Four of them are determined by how much a(r) increases over each stage based on 
the cosmological considerations. As in [29, 30, 31], we take the following specifications: ^1 = = 300 
for the reheating stage, (s = = 10^"^ for the radiation stage, C2 = = 3454^^^ for the matter 
stage, and = a(T°J) ~ iu^)^^^ present accelerating stage. The remaining time instance is fixed 

by an overall normalization 

|ro - Ta\ = 1. (8) 

Notice that this convention of normalization is different from that of Ref. [19]. There are also 12 constants 
in the expressions of a(r), among which /3, (3s and 7 are imposed as the model parameters, for the inflation, 
the reheating, and the acceleration, respectively. Based on the definition of the expansion rate Hq = ^\to 
of the present universe , one has Ih = 'j/Hq. Making use of the continuity conditions of a(r) and of 
a(T)' at the four given joining points ri, Tg, T2 and te, all parameters are fixed as the following: 

Ta-TE = Ce) 
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T2-Te = -C2 'Cs' 

7 

1 _i i 

7 

r,-rp = -(l + /3,)C^C2"^c|, 
7 

n - Tp = i(i + /3.)Ci'^C'C2""c|, 

ri = ^(l + /3)Cr^C'C2"'c|, 

Te^ = 2.75 X 10-" (r2 - Te) - Te, (9) 



— ^ 7 SB ) 

Oe = 7 C2 ^ Cb , 

^fs-l -{l+l±ik) 



h = iH7''-''\i+p\-^'^''ki'^''Ck^c~E ' ■ (10) 

In the expanding Universe, the physical wavelength A is related to the co-moving wave number k by 

and the wave number corresponding to the present Hubble radius is 

J-Z-no 

There is another important wave number involved 

_ 27Ta{TE) ko , . 

kE=-rjj— = ——, (13) 

whose corresponding wavelength at the time te is the Hubble radius I/Hq. In the present universe the 
physical frequency u corresponding to a wave number k is given by 

= — = — - = k. (14) 

A 27ra(ro) 27r7 



3. Analytical solution of RGW 
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In the presence of the gravitational waves, the perturbed metric is 

ds'^ = a'^{T)[-dT'^ + {Sij + hi-j)dx^dx^], (15) 
where the tensorial perturbation hij is a 3 x 3 matrix and is taken to be transverse and traceless 

h\ = 0, h,,,, = 0. (16) 

The wave equation of RGW is 

a,(^9^/ty(T,x)) =0, (17) 

where g is the determinant of the 4-dimensional metric y^^. One can decompose hij into the Fourier 
modes as 

h^J{r,^) = E/ ' (18) 



where k is co-moving wave vector, a denotes polarization states +, x, and e^y is the polarization tensor 
Since RGW is isotropic, Eq. (17) reduces to 



%"(r) + 2^hk'iT) + ehkir) = 0. (19) 

a(r) 

where the co-moving wave number k is the Euclidean length of k. For each polarization, x, +, the 
equation of h^^^ is the same, so the super index {a) can be dropped from h^^\ Since for all the stages of 
cosmic expansion the scale factor can be approximately described by a power-law form 

a(T) oc r", (20) 

the solution to Eq. (19) is a linear combination of Bessel function J^, and Neumann function N^, 

hk{T) = x^^/^^-''[aiJ^^^i/2){kT) + a2N^-ii/2){kT)l (21) 

where the constants oi and a2 can be determined by the continuity of and of h'j^ at the joining points 
7"i,rs,r2 and te- However, as mentioned above, during t^w < t < T2, Eq. (19) will be modified and its 
solution will be given later. 

The inflationary stage has the solution 

%(r) = Aol^M~^^'^^'^[MJii/2)+p{p) + A2J-((i/2)+/3)(p)], -oo < r < n (22) 
where p = A;r and 

M = V^/l^^'^^^'' ^^ = ^\/f^"^'^^^'' (23) 

cos pvr V 2 cos pvr V 2 

are taken from [33], and thus the so-called adiabatic vacuum is achieved: liuik^oo hk(T) oc e^*'^"^ in 

the high frequency limit [34]. The constant Aq determined by the initial amplitude of the spectrum is 
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independent of k. For kr ^ 1, hk^r) oc J(i/2)+i3ix) oc As seen below, this choice will lead to 

the required scale-invariant initial spectrunn in Eq. (41). 
The reheating stage has the solution 

%(r) = t-(V2)-/3. [Si J(i/2)+;3, (A; t) + B2N^y2)+p, {k t)], n < r < r, (24) 

where t = T — Tp and 

Bi = -^7rtp)+^nfeiV(3/2)m(^«i)^fc(ri) + iV(i/2)+/3.(A;ti)4(ri)], (25) 
B2 = ^7rif/')+^^[fcJ(3/2)+^,(A:ii)/ifc(Ti) + J(i/2)+/3.(A;ti)4(Ti)], (26) 

with ti = n — Tp, and /ife(ri) and /i]fc(Ti) are the corresponding values from the precedent inflation stage. 

The radiation-dominant stage needs to be divided into two parts. The first part of this stage is before 
the EPT when < r < Tew, the PMF has not been produced yet, so the wave equation of RGW is still 
homogeneous with the solution 

hk{r) = [CiJi/2{ky) + C2Ny^{ky)] , Ts<T<Tew (27) 

where y = t — Tg and 

Ci = -l^ys^''[kNs/2{kys)hk{Ts) + Ny2{kys)K{Ts)i (28) 

C2 = \'Kyl'^[kJz/2{kys)hk{Ts) + Ji/2{kys)h'k{Ts)l (29) 

with ys = Ts — Te, and hk{Ts) and h'^{Ts) are from the reheating stage. The second part is from the EPT 
up to the matter-dominant stage when Tew < t < T2. During this period the wave equation of RGW is 
modified as 

h'lAT) + 2^^h\AT) + k^hij{T) = 16TrGa''Uij, (30) 

a[T) ■' 

where Iljj is an anisotropic term contributed by the PMF. Note that the righthand side of Eq. (30) is 
different from the Eq. (15) in Ref. [19], this is because we choose a different normalization of scale factor 
|ro - Ta\ = 1. As in [19], setting n(A;,r) = (^)2/(A;, r)n(fc) (where oq is today's scale factor and Ii{k) 
is a time independent random variable with power spectrum (|n(fc)|^) = 1), Eq. (30) can be reduced to 

hl{T) + 2^/i',(^) + k^hk{T) = 167rGalf{k,T). (31) 

For the PMF produced by EPT, the source is a white noise independent of k [19], and the key function 
/ can be written as (See Appendix A for detailed derivations) 

/(^)-%;^M.)-<-. (32) 

Bx is the amplitude of PMF, A is the length scale on which cosmic magnetic fields have been observed. 
For the PMF produced by EPT, the spectral index n = 2, and kc = i/iiin with rjin — 4: x lO^sec is the 
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corresponding cutoff scale [19]. By using the theorem of Wronskian, Eq. (31) can be analytically solved 
(See Appendix B for details), and then the effects of PMF can be included into our calculation of RGW. 
The nnatter-dominant stage has the solution 



%(r) = 



, r2<T <TE 



where q = T — Tm and 



Di = -\'r^<iT[kN^/2{k<.2)hk{T2) + N^/2{k<.2)h'k{T2)], 
D2 = \'^4!'^[kJ^/2{k(.2)hk{T2) + Jz/2{k<;2)h'k{T2)\, 



(33) 

(34) 
(35) 



with <j2 = T2 — Tm. In the expressions of Di and D2, the mode functions hk{T2) and h'j^{T2) are also from 
the precedent stage. 

The accelerating stage has the solution 



where s = r — r„, and 



,(l/2)+7 



£^1 J_(l/2)-^(A: S) + £^2^-(l/2)-7(fc s) , Te <T <To 



Ei = —--K s 



(l/2)-7 



[^^(l/2)-7(^ SE)hk{TE) + N_^y2)-^{kSE)K{TE) 
-^2 = \'!^ S^^^^^ [kJ(j^l2)-y{k SE)hk{TE) + J-{l/2)-^{k SE)h'k{TE)], 



(36) 

(37) 
(38) 



with se = Te — Ta- So far, the explicit solution of /ife(r) has been obtained for all the expansion stages, 
from Eq. (22) through Eq. (36). 



4. Theoretical Predictions of RGW 

There are three different spectra that are commonly used to describe the stochastic GW background: 

the tensor power spectrum A^(A;, r), the gravitational wave spectrum /i(/c,r), and the energy density 
spectrum Q.gyj{k,T). In the early universe, the RGW is usually characterized by A^(/c,r); on the other 
hand, the present-day RGW is usually characterized by h{k,T) or by Q.gyj{k,T). In particular, as pointed 
by Grishchuk [25], the spectrum of RGW h{k,T) can be directly compared with the direct measurements 
of Laser interferometer GW detection, and is defined by the following equation [25]: 

POO fJU 

h\k,T)^^m^{^,T)h,{^,T)\0), (39) 
where the right-hand side is the expectation value of the h^^hij. Calculation yields the spectrum at present 

Mfc,ro) = -fc3/>fe(To)|. (40) 
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One of the most important properties of the inflation is that the initial spectrum of RGW at the time Tj 
of the horizon-crossing during the inflation is nearly scale-invariant [25]: 

/ k \ '2+13 

h{k,n) = A(-) , (41) 



where constant A is directly proportional to Aq in Eq. (22), and will be fixed by the observed CMB 
anisotropies in practice. Since the observed CMB anisotropies [35] are AT/T ~ 0.37 x 10~^ at I ~ 10, 
as in Refs. [30, 31], we take the normalization of the spectrum 

h{kE, To) = 0.37 X IQ-^r^ , (42) 

where ks is defined in Eq. (13), and its physical frequency being ue = A;f;/27ra(ro) = Hq/{1 + ze) ~ 
1 X 10~^^ Hz. The tensor/scalar ratio r can be related to the slow roll parameter e in the scalar inflationary 
model as r = 16e [26]. However, the value of r is model-dependent, and frequency-dependent [36, 37]. 
This has long been known to be a notoriously thorny issue [38]. In our treatment, for simplicity, r is 
only taken as a constant parameter for normalization of RGW. The WMAP5 data combined with Baryon 
Acoustic Oscillations and Type la supernovae give r < 0.20 (95% CL) [27]. For concreteness, we take 
r = 0.20. 

Next, we shall focus on how the theoretical prediction of RGW is related to the GW experiment. As 
pointed by [39], 

(0|/i*^(x,ro)/ii,-(x,ro)|0)=4/ uSh{v) — , (43) 

where Sh{u) is the one-sided spectral density of RGW. Base on Eqs. (39), (43), and using Eq. (14) to 
change dk/k into dv/v, one gets 

The noise level of the detector is measured by the strain sensitivity /i„, defined as 

K = ^JSniiy), (45) 

where Sn{i') is the one-sided spectral density of the detector's noise. A stochastic GW background will 
be observable at a frequency v if [39] 

Sh{v) > Sn{v)/F, (46) 

where F is the angular efficiency factor determined by the geometry of the detector. For interferometers, 
F = 2/5 is always satisfied [39]. Therefore, to compare the theoretical model of RGW with the latest 
LIGO data, one needs to calculate 

^(^'^o) (47) 



VlOz/ 

In principle, since all the RGW solutions of the 6 cosmic stages are analytic, the analytic expression of 
/i(i/, To) can be written down. However, since this expression is too complex, we only give the numerical 
results in this work. 



9 




-20 -15 -10 -5 5 10 



log,„ V [Hz] 

Figure 1: The spectrum h^v, tq) of RGW for various amplitudes Bx of PMF. PMF will not markedly 
affect h{v,To), unless its amplitude is stronger than Bx > 10~^ Gauss. 

5. Results and Discussions 

In the following, by comparing our theoretical results with LIGO S5 data, we will give the constraints 
on the amplitude of PMF. 

In Figure 1, we plot the spectrum h{i', tq) of RGW for various amplitudes Bx of PMF, where the scale 
of galactic cluster A = 1 Mpc is adopted. PMF yields a level spectrum of RGW in the range (10^^ — 10'^) 
Hz, which covers the bands of operation of LIGO. For a larger amplitude Bx, this level spectrum of RGW 
will also be higher. It is seen that PMF will not markedly affect tq), unless its amplitude is stronger 
than Bx > lO^^ Gauss. 

Figure 2 is a comparison of the theoretical predictions of RGW with the latest LIGO S5 data [24]. 
Again, a scale of galactic cluster A = 1 Mpc is adopted in this comparison. The three fluctuating curves 
are the RGW's theoretical prediction of Bx = 10^^ Gauss, Bx = 10^^ Gauss, and Bx = 10^^ Gauss, 
respectively. The dashed line is the latest lower limit given by the LIGO S5 run [24]. Here the vertical axis 
is the strain sensitivity, given by Eq. (47). We find that for the scale of galactic cluster A = 1 Mpc, the 
amplitude of PMF has to be weaker than -Ba < 4 x 10^'' Gauss. 

Now, let us turn to the case of PMF with a scale of supercluster. Figure 3 is a comparison of the 
theoretical predictions of RGW with the latest LIGO S5 data, where the scale of supercluster A = 100 
Mpc is adopted. The three fluctuating curves are the RGW's theoretical prediction of Bx = 10^^ Gauss, 
Bx = W^^^ Gauss, and Bx = 10^^^ Gauss, respectively. The dashed line is the latest lower limit given 
by the LIGO S5 run [24]. It is found that for the case of A = 100 Mpc, the amplitude of PMF has to be 
weaker than -Ba < 9 x 10^^^ Gauss. 

In this paper, we present new PMF limits by using LIGO S5 data. Compared with the constraints from 
the nucleosynthesis limit [19], our results are much closer to that from the WMAP5 -B-mode polarization 
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Figure 2: Comparison of the theoretical predictions of RGW with the latest LIGO S5 data, where 
the scale of galactic cluster A = 1 Mpc is adopted in this comparison. 




Figure 3: Comparison of the theoretical predictions of RGW with the latest LIGO S5 data, where 
the scale of supercluster A = 100 Mpc is adopted in this comparison. 
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limit [20]. Although our constraints are not as strong as those from CMB observations, future laser 
interferometer GW detectors, such as advanced LIGO [40], BBO [41] and DECIGO [42], will be able to 
put much tighter constraints on the PMF. In this manner, current and future GW detectors have potential 
to make interesting contributions to the study of the PMF. 
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Appendix A: Detailed derivations of the key function / 

As in Ref. [19], we shall only consider the magnetic fields on sufficiently large scales, and ignore the 
impacts come from small scales. We model today's magnetic field Bo(x) as a statistically homogeneous 
and isotropic random field. The transversal nature of magnetic field B then leads to 

(5^(k)S*^(q)) = d^{k - ci){d'^ - Uy)B'^{k), (48) 



where k and q are co-moving wave vectors, k = 'k/k, and k = ^Jj2iik^)'^ is the Euclidean length of k. 
We use the following Fourier transform conventions 

B^{k) = ^exp(ix • k)5^(x), S^(x) = ^ f d^kexp{-ix ■ k)B^{k), (49) 
k 

where V = J d^r cxp(—j^) is the normalization volume (A is the length scale of cosmic magnetic fields). 
On sufficiently large scales, ^^(x) = (^^)^S^(x), 5,(x) = (^|i)So,(x), and 5*(x) = (^^y^5^(x). If B 
is generated by a causal mechanism, it is uncorrelated on super horizon scales, 

(B*(x,r)5^ (x',r)) = for |x - x'| > 2r. (50) 

Notice that the universe is in a stage of standard Friedman expansion when the PMF is produced, so that 
the co-moving causal horizon size is about r. According to Eq. (50), (5'(x, r)B-'(x', r)) is a function 
with compact support and hence its Fourier transform is analytic, i.e. 

{B\k)B*^{k)) = {6'^ - k'k^)B^{k) (51) 

is analytic in k. The Maxwell stress tensor T^^ of a magnetic field in real space is given by 

r^-(x,r) = ^\b%^,t)B^{^,t) y^{yi,T)Bi{:K,T)B\^,T)\. (52) 
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In Fourier space, making use of Eq. (49) and the scaling of the magnetic field with time, we have 

r^(k,r) = i-(^)'^ [5^(q)5^(k- q) - ^Bi{ci)B\k - ^1)6'^. (53) 

n*-'(k, r) is the transverse traceless component of r'^(k, r) and sources the RGW, defined as 

ff^(k, r) = {P^Pl - ip^^P„6)T"''(k, r), (54) 

where = Sij — kikj is the projector onto the component of a vector transverse to k, and P^P^ projects 
onto the transverse component of a tensor. Here we give the details of the calculation of its correlation 
function (n*-^(k, r)n*'™(k', r)), which will be used to compute the induced GW. Defining the projector 

'P^ = PaPi - iP^'Pab, we have 

(ff^(k, r)n*'"*(k', r)) = V'\,,V^1i{T''\k, T)T*'^{k', r)). (55) 
Note that up to a trace, which anyway vanishes in the projection (55), T"''(k, r) is just given by 

A^\k, r) ^ (^) ' ^ B-{ci)B\k - q). (56) 



47r V a 
Therefore, we can get 

(ff^(k,r)n*'™(k',r)) = P\,P^™,(A«''(k,T)A-'i(k',T)). (57) 

Assuming the random magnetic field be Gaussian, the Wick's theorem for Gaussian fields can be applied 
here. Making use of Eq. (48), Eq. (56), and the reality condition B*°'(k) = B"'{—k), products of four 
magnetic fields can be reduced as 

(P«(q)P''(k-q)P*'=(p)S*<^(k'-p)) 

= Wk(<5"' - rq')B\q) ■ <5-p,k'-p(<^^' - rq'')B\-p) 

+5q,p(,5- - rr)B\q) ■ 4-q,k'-p('5'"' - k^qk^q^)B\\ci - k\) 

+'5q,k'-p(5"' - rq'')B\q) ■ 5k-q,p(<5'^ - k'^q'k'^q'')B\\q - k|). (58) 

The first term only contributes an unimportant constant and then can be disregarded. Making use of Eqs. 
(56), (57), (58), and setting Vabcd = 'PijabPcd = '^aUj'P]^, one can obtain 

(n,,(k,T)n«^(k',T)) = (-)X,ed(A"^k,r)A-'='(k',r)) 

ao 

= E B\q)B\\k - q|)(l + 27^ + 7^/3^), (59) 

q 

where 7 = k • q, and /3 = k • k — q. As in Ref. [19], setting 



(n,,(k,r)n*^^(k',r)) = 4(^)V'(A;)<5k,k', (60) 
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one can get 

^'^^^ = ^64^^ ^ ^'(5)^'(|k - q|)(l + 27' + 7'/3'). (61) 
q 

In a coordinate system where k is parallel to the z-axis, 11^ has the form 

/ TT TT n \ 

Y u u u y 

where + and x denote the polarization states of GW. Base on Eq. (60), the statistical isotropy gives 

m+\') = {\^.\') = i-rf- (62) 

a 

To continue, we have to specify B'^{k). Base on Eq. (49), one gets 

(i?o(k)i?o(x + r)) = 1 / d3^^^e-(''+''')-''e-'^'-"-i?(k)S(k'), (63) 

k k' 



n+ 


Hx 







-n+ 















and the normalization is 



Bl = ^J ciV(Bo(k)5o(x + r))e-^. (64) 

As in Ref. [19], we assume that B'^{k) can be approximated by a simple power law, i.e. B^{k) = Ck^ (C 
is a constant). Base on Eq. (64), a tedious but straight forward computation gives 

_ Bl{X/V2r 

^-27r5/2r(i^)' ^^^^ 



and thus 



,2,.^ _ Bl{X/V2r 



B\k) = , 7 ' (66) 
^ ' 27r5/2r(2^) ^ ' 



Making use of Eqs. (61) and (66), one gets 

/'(^) = slniff^^) / d\{in\^^ - q|)"(l + 27^ + (67) 

Now let us focus on the PMF produced by EPT. As mentioned above, for this case, and the source is a 
white noise independent of k. Utilizing the integration result of Ref. [19], finally we get 

, _ Bi{X/V-2r^^ 

where the spectral index n = 2, and the corresponding cutoff scale is kc = i/rjin with rjin — 4: x lO^sec 
[19]. This leads to the expression of key function / we need in Eq. (31). Note that the expression of / 
has subtle difference with that in [19], this is because we use different conventions of dimension. 

Appendix B: Exact Analytic Solution of Differential Equation (31) 
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As mentioned above, to obtain the solution of RGW during the period T(,w < t < T2, one must solve 
the differential equation (31). Follow [19], we can write the mode function as 

hkir) = hk{Tew)x{u), (69) 

where hk{Tew) is given by Eq. (27) evaluated at Tgw, and x(u) satisfies the following differential equation 

,",„) + 4» + xW = ^^. (70) 

with u= kr, and s{k,T) = 167rGaQ/(A;, r). The homogeneous solution of Eq.(70) is 

, , cos(u) sin(u) 

u u 

where the coefficients 

Wi = Uew COs{Uew) + {WUew - 1) Sm{Uew), (72) 
W2 = Uew Sm{Uew) - {WUew - 1) COs(Ue,„), (73) 

are fixed by the continuity condition of % and h'f^ at the time instance r = Tew, with Uew = krew, and 

^ _ BiJ3/2[kiTew - Te)] + ^2-/V3/2 [^('T'e^) - Te)] 
BlJi/2[k{Tew - Te)] + B2Ni/2[k{Tew - Te)] ' 

Using the theorem of Wronskian, one can get a special solution of Eq.(70) 

[/(sin li cos [/ - sin [/ cos «) rr^JTr rr,r\ 

Therefore, the general solution of Eq.(70) can be written as 

x{u) = Xi{u)+X2iu), (76) 
and thus the solution of /^^(t) at the period Tew < r < T2 is in hand. 
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